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Abstract 

The fractional Fourier transform (FRFT) naturally exists in the strongly nonlocal nonlinear 
(SNN) media and the propagation of optical beams in SNN media can be simply regarded as a self- 
induced FRFT. Through FRFT technic the evolution of fields in SNN media can be conveniently 
dealt with and an arbitrary square-integrable input field presents generally as a revivable higher 
order spatial soliton which reconstructs its profile periodically after every 4 times of Fourier trans- 
forms. The self-induced FRFT would illuminate the prospect of the SNN media in new applications 
such as continuously tunable nonlinearity- induced FRFT devices. 
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Since Fourier suggested the usage of Fourier analysis to solve the heat conduction problem 
in 1807, the Fourier transform (FT) has been applied widely in many branches of science [l| ■ 
In the field of optics, the FT is one of the most important and basic tools in dealing with 
physical optics and optical information processing l|. In fact, the term Fourier optics is 
often used synonymously with optical information process. The fractional Fourier transform 
(FRFT), which is an extension of the FT {2J], has been introduced to optics since 1993 when 
Mendlovic and Ozaktas find this operator can be optically performed by the quadratic 
graded-index (GRIN) media [3]. In respect that the FRFT can show the characteristics 
of the signal changing continuously from the spatial domain to the spectrum domain, it 
interests optics scientists and engineers and plays an important role in man y op tics fields 



llQl: 



12], such as diffraction [7J], transmission 8|, |9j, imaging lOj, information 



processing LUJ, LL^IU etc. 

On the other hand, since the first observation of nonlinear optical phenomena by Franken 
et al in 1961 jl3j |. one year after the invention of the laser, the nonlinear optics has been 
a rapidly expanded active field and widely influenced other fields. A special branch of the 
nonlinear optics, the field of optical solitons, have grown enormously in the past decades. In 
particular, the soliton in strongly nonlocal nonlinear (SNN) media which supports (2+1)- 

3een widely investigated in the 



dimension solitons, have attracted extensive interest and 
past few years 15, 16, 17, Q, Q, 21, 22, 23, 24, 25, 26J 



271 ]. since Snyder and Mitchell 



simplified the nonlocal nonlinear Schrodinger equation (NNLSE) to a linear model (called 
Snyder-Mitchell mode (SMM) in our papers) in the SNN case and found an exact Gaussian- 
shaped "accessible soliton" [14] . And the "accessible soliton" has shown its interest in po- 
tential applications such as photonic switching and logic gating 27]. 

The Fourier optics and the nonlinear optics might seem to be independent of each other 
in that the FT and the FRFT are linear transforms and well known as technic provided 
to solve problems in linear system. But we note that the SNN media, the propagation 



equation in which can be mathematically simplified to the linear SMM 14j, would provide 
an opportunity to intersect them with each other. 

In this letter, the FRFT induced by the input field itself through the SNN effect is 
introduced. It is shown that the FRFT naturally exists in the SNN media. The order of the 
self-induced FRFT in the SNN media can be steered with the input power in addition to the 
propagation distance z, quite different from the traditional linear (intensity-independent) 
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FRFT devices such as the lens series and the quadratic GRIN media. Due to the self- 
induced FRFT, propagation of an arbitrary square-integrable input field presents generally 
as a revivable higher order spatial soliton (RHOSS) which is similar to the higher order 
(1+1)D temporal soliton in fiber. The build-block-like solitons, build-block-like breathers, 
and their interaction, can be regarded as special cases of the RHOSS and can be conveniently 
investigated with the FRFT technic. 

The propagation of beams in nonlocal nonlinear media can be phenomenologically de- 
scribed by the NNLSE: 2ikriod z $> + n A^$ + 2£; 2 An$ = 0, where k represents the wave 
number in the media with the linear part of the refractive index no when the nonlinear 
perturbation of refractive index An equals zero, An = n 2 J R(f — r a )|$| 2 d 2 r a (n 2 is the 
nonlinear index coefficient and R is the normalized symmetric real spatial response function 
of the media). In the case of SNN media we need only keep the first two terms of the 
expansion of An and the NNLSE is simplified to the SMM 14| 

2ikd z <$> + A ± $ - A; 2 7 2 P r 2 $ = 0, (1) 

where 7 is a material constant, Pq = J |$| 2 d 2 r is the input power. 

We assume $(r, z) = a(r) exp(— i/3z) to seek the stationary solutions of Eq. ([I]). Substi- 
tuting this expression into Eq. gives 2/3ka = k 2/ y 2 Por 2 a — A±a, the eigen solutions of 
which are Hermite Gaussian soliton family 

a® = HG® n = c^H m (—)H n {^-)e~A (m + n = p) (2) 

w c w c 

in cartesian coordinate [26], Laguerre Gaussian soliton family 



. r 2 
2w'i 



lg c(p) = c (p),j L y L i ( JL \ cos no) e 
a^={ l ' q KwJ qKwJ 1 _ 2 (2q + l=p) (3) 

LG^ = c^\0L^)sm(W)e-^ 



in circular cylindrical coordinate 26] , or Ince Gaussian soliton family in elliptical coordinate 



23|, 



241 ] . In Eqs. ([2]) and (j3J), p = 0, 1, 2... is the order of the solution with which the soliton 
eigenvalue can be straightforwardly obtained: 

/?(p) = (p + l)/3b, (4) 

is the normalized coefficient ensures J |a^| 2 d 2 r = Pq, w c = (A; 2 7 2 Po)~ 1 ^ 4 is a gener- 
alized beam width, H m (x)a,nd L\ Jr) respectively represents the Hermite and the Laguerre 
polynomials, j3 = y/Poj = 1/kw 2 . 



Mathematically, the profile of the eigen soliton solution, aW, is simultaneously the eigen 
function of the FRFT |a, hi, yj 

FMn)m = 7' <( °:'" e»p[ r J±-] /«p[ * M MiW , 

27ru;^sina 2^ tan a J 2^ tan a u^sma 

i.e., is the eigen solution of the equation F a {a^(fi)} = a^ p \f2)e~ ipa with the eigen 
value e~ wa , where the order of the FRFT is 



a = (3 z = a/ Pq^z. (5) 

Therefore the field of the eigen soliton at z is connected with that at the entrance plane 
through the FRFT: ¥P\r 2 ,z) = F a {^{f 1} 0)} x exp(— ia). 

An arbitrary square-integrable input field can be expressed as a linear superposition of an 
arbitrary one of the three families of eigen soliton solutions: $(f 1; 0) = c p $ (p) (n, 0). 

According to the linearity of Eq. ([I]) and the FRFT, the propagation in the SNN media is 
presented as the FRFT on the input field 

$(f 2 , z) = 0)} x exp(-ia). (6) 

In the special case a = ir/2, the propagated field at z is deduced to the well-known FT of 
the input field, reads 



$(r 2 , z) = — — / $(f x , 0) exp[ L T ^]d 2 f 1 . (7) 



Equations (E])- (JTJ) indicate that the SNN media naturally performs the FRFT and FT on 
the input field. The physical origination of this effect is as follows: When a beam is input 
into a SNN media, it would induce a quadratic GRIN channel in the medium through the 
nonlocality. The propagation in the channel then performs the FRFT and FT, just as in the 
traditional quadratic GRIN medial . Because the grads of the refractive index distribution 
can be steered with the power, the order of the self-induced FRFT in the SNN meida can 
be steered with the input power in addition to the propagation distance z, quite different 
from the traditional linear FRFT devices. 

The simulation based on the NNLSE shows that: the self-induced FRFT becomes more 
and more distinct with the increase of the nonlocality or the input power (Fig. (He)). In 
fact, the increase of the power is tantamount to the increase of the nonlocality, in that it 
decreases w c and the average beam width in propagation. On the other hand, the smoother 
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FIG. 1: (Color online) Porpagation dynamics of the super Gaussian field exp 
media with Gaussian response function R(r) = l/(2irw R ) exp[— r 2 /(2w^)] 
numerical simulation of the NNLSE. (a) Evolution of the profile in propagation, (b), (c) Intensity 
distribution at f3$z = tt/2 and it respectively. The nonlocality degree T = wb/wr = 1/10, wb is 
the second-order moment width of the beam at the entrance plane, P\ is the critical power for the 
eigen soliton with generalized width w\. 

the shape of the input field is, the less the nonlocality is required to support the FRFT, 
because it contains less higher-order eigen soliton solutions in the superposition. Generally, 
the profile of the eigen soliton holds when V < 1/10. When the order of a constituent eigen 
soliton in the superposition is high enough so that V > 1/10, the profile would be distorted 
in propagation. Therefore the SNN media act as a low-pass spatial frequency filter. 

Based on the self-induced FRFT, it is convenient to investigate the propagation in SNN 
media from the angle of Fourier optics. We can predict the field would present a periodical 
evolution with the period Az = 27r/y / ^)7 (corresponding to Aa = 2n) (Figs. OB). As 
a result of cascade FT, the pattern evolves to the inversion of the input pattern at z — 
(2n + l)7r/v^7 and revive to the input pattern at z — 2mr/y/Po'y (Fig. [2]) (we call those 
cross sections the revived planes or imaging planes). At z = (n+ l/2)ix / y/T^ff (we call these 
cross sections the Fourier planes) , the patterns are the FT spectrum of the input field or the 
inversion. From Eq. (jTJ) the spatial frequency k r = f^/w 2 , thus at the Fourier planes the 



5 



p Cl z=0 7/4 a/2 3n/4 a 




5^/4 3*/2 7*/4 2ji 9^4 




FIG. 2: Propagation dynamics of Gaussian beams (exp[— (x/\/2wi) 2 ]) truncated by triangular 
(rows 1 and 3) and square (rows 2 and 4) super Gaussian diaphragm in SNN media with Gaussian 
response function based on numerical simulation of the NNLSE. Pq = 4Pi. T = 1/20. 

beam width w(z) oc l/y^Po, which obeys the scaling rule of FT (Fig. [0(b)). Because of the 
periodic revivable evolution, which is similar to the higher order temporal soliton in fiber, 
we call this type of propagation the RHOSS (The RHOSS should be distinguished from 
the traditionally mentioned "higher order spatial soliton" which refers to the stationarily 
propagated multipole spatial soliton). There is an interesting difference between the higher 
order temporal soliton in fiber and the RHOSS: the existence of the higher order temporal 
soliton in fiber requires much more power than that required for the fundamental soliton, 
whereas to support the RHOSS, the power can be lower than the critical power of the 
stationarily propagated fundamental soliton. 

According to properties of the FRFT, there are three special cases of the RHOSS (Fig. 

H: 

1) Build-block-like soliton. As shown in Eqs. (|4|), during propagation, every degenerate 
eigen soliton with the order p has the same propagation-induced phase (3^z (or in other 
words, has the same FRFT eigen value exp(— ipa)). Therefore, when i) the input field is the 
linear superposition of the eigen soliton solutions with the same order p, beam center, and 
generalized width wi; and ii) the power is the critical power Pi = 1 / \k 2 w\'-) 2 ) which supports 
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the eigen solitons with the generalized width w±, the field would propagate stationarily, i.e., 
the soliton occurs (rows 1-2 in Fig. [3]). Because the field of this type can be freely composed 
of the eigen soliton solutions with the same order, we call it build-block-like soliton. 

2) Build-block-like breather. When all are the same as the build-block-like soliton except 
the input power P deviate from the critical power P x , the FRFT keeps the shape of the 
input field and periodically varies the width with the period Act = n. Subsequently the 
evolution is presented as breather with the period Az = Tr/v^rT (row 3 in Fig. [3]), and 
the change of the build-block-like breather's width is the same as the prediction for the 
Gaussian breather in Ref. [uj]. At z = (n + l/2)7r/v/^)7, the FRFT is deduced to the FT 
and w(z) oc l/y/P^w\. 

3) Multi soliton interaction. In this case, the traditional technic might require much effort 
in mathematical treatment because orbits of the interacting solitons are influenced by each 
other. But by introducing the FRFT, the evolution of the orbits is presented simply as 
a shift in the FRFT reads <&(r2,z) = -F a {$(r*i — ro, 0)} exp[— ia], where r*o represents the 
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FIG. 3: (Color online) Propagation dynamics of build-block-like soliton (rows 1-2), build-block- 
like breather (row 3), and multisoliton interaction (row 4) in SNN media with Gaussian response 
function based on numerical simulation of the NNLSE. T = 1/10, Pq = 2P\ for row 3 and Pq = P\ 
for others. The input fields are respectively (HGf\ - FG^/10) (row 1), (HG^l + HGfl) (row 2), 



(LG^f +iLG { 8 y> -HG^/125) (row 3) and [HG$(x+toi,y+^i)+HG%{ ) (x-3wi,y-&wi)](Tcm 
4). The general width of all constituent beams are w\ at the entrance plane. 
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FIG. 4: The sketch of the experiment setup. 




FIG. 5: Experimental results demonstrating self-induced FRFT of Gaussian beams truncated by 
the triangular ((a)-(f)) and the square ((g)-(l)) diaphragm at different input power, (a) and (g) 
are the input fields, (b)-(f) and (h)-(l) are the output patterns. 



initial deviation of the beam center of the interacting soliton from the mass center. Under 
the vertical incidence condition (row 4 in Fig. [3]), the solitons intersect each other at z — 
(n + 1/2)^/^/^7, evolve to the inversion $(r*o — r, 0) at z = (2n + l)7r/v^7, an d recur to 
the input field $(f — r , 0) at z = Iwn j y/P^ . 

To verify the prediction about the self-induced FRFT (or in other words, the RHOSS), we 
carried out the experiment in a columned lead glass. The experimental setup is illustrated 
in Fig. HI The beam from a Verdi laser is focused by the collimation lens. A diaphragm is 
placed at the focus and the real image is produced at the entrance plane of the lead glass 
by the confocal lenses pair. When the input power is adjusted, the intensity distribution at 
the entrance and exit plane are monitored by imaging the beams onto a CCD camera. The 
lead glass is 59.8 mm in length and 15.1 mm in diameter. The beam at the entrance of lead 
glass is 85 /mi x 80 /im in size (for square diaphragm) or 106 /im in diameter of circumcircle 
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(for triangle diaphragm). 

The experimental results are shown in Fig. [5fa)-(l): by changing the input power, the w c 
changes and the patterns similar to the FRFT spectrums with the orders a = 0, 7r/4, 7r/2, 
37r/4, and 7r in Fig. [2] are recorded. In case of the triangle (square) diaphragm, the output 
pattern recur to the input pattern when P = 551 mW (Fig. E^b)) (P — 590 mW (Fig. 
Efh))) and evolves to the inversion when P = 823 mW (Fig. [5]^f)) (Po = 810 mW (Fig. 
EJ4))). Because the higher-frequencies are filtrated in propagation, the output patterns are 
smoother than the input ones. In Fig. [5](z), the variation of the FRFT order a with the 
square of the input power, i.e. \^P~o, is illustrated. The linear fit shows that the FRFT order 
a is directly proportional to \/Po, as predicted in Eq. ([5]). 

In summary, the self-induced FRFT in SNN media made the nonlinear optics and the 
Fourier optics intersect with each other. The introducing of the FRFT technic would release 
one from complicated mathematical calculation in propagation problems such as soliton 
solutions in SNN media. The RHOSS, including the build-block-like solitons, build-block- 
like breathers, and their interaction, would greatly enrich the nonlocal soliton family. The 
fact that the order of the self-induced FRFT in SNN media is related not only to the 
propagation distance but also to the input power, quite different from that in the traditional 
linear devices, would illuminate the prospect of new applications of the SNN media such as 
developing power-controlled continuously tunable FRFT devices. 
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10674050 and 10804033), the Program for Innovative Research Team of the Higher Education 
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